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Afterglow, or long-lived emission, has now been detected from about a dozen 
well-positioned gamma-ray bursts. Distance determinations made by mea- 
suring optical emission lines from the host galaxy, or absorption lines in the 
afterglow spectrum, place the burst sources at significant cosmological dis- 
tances, with redshifts ranging from ~1— 3. The energy required to produce the 
bright gamma-ray fiashes is enormous: up to ~ 10^"^ erg or 10% of the rest mass 
energy of a neutron star, if the emission is isotropic. Here we present the dis- 
covery of the optical afterglow and the redshift of GRB 990123, the brightest 
well-localized GRB to date. With our measured redshift of > 1.6, the inferred 
isotropic energy release exceeds the rest mass of a neutron star thereby chal- 
lenging current theoretical models for the origin of GRBs. We argue that the 
optical and IR afterglow measurements reported here may provide the first 
observational evidence of beaming in a GRB, thereby reducing the required 
energetics to a level where stellar death models are still tenable. 



Almost thirty years after their discovery^, the distance scale to gamma-ray bursts 
(GRBs) was finally unambiguously determined just two years ago^. Thanks to precise 
and prompt GRB localizations by the Italian-Dutch satellite, BeppoSAX^ and the All 
Sky Monitor^ aboard the Rossi X-ray Timing Explorer, afterglow emission spanning the 
wavelength range from X-ray to radio has now been detected in about a dozen events^'^'^. 
This has provided redshift measurements for four bursts: one through the detection of 
optical absorption lines^ and three by identifying emission features in the host galaxies^'^'^'^. 
It is now assumed that a substantial fraction of all GRBs occur at significant cosmological 
distances {z ~ 1-3). 

The large inferred distances imply staggering energetics. In particular, our measured 
redshift of z = 3.14 for GRB 971214 (ref. 8) led us to an inferred isotropic energy loss of 
3 X 10^'^ erg, or O.IMqc^. At that time, this result was seen to challenge most proposed 
GRB models. The energy estimate can be reduced by invoking non-spherical emitting 
surface geometry, e.g. jets. However, to date, there is no firm observational evidence for 
jets in GRBs. Indeed, on statistical grounds, Grindlay^^ and Greiner et al.^^ constrain the 
"beaming factor" (the reduction in inferred energy release relative to that for the spherical 
case) to be less than 0.1. 

BeppoSAX ushered in 1999 with the discovery of GRB 990123^^, the brightest GRB 
seen by BeppoSAX to date. It is in the top 0.3% of all bursts, if ranked by the ob- 
served fiuence^^. If the inferred energetics of GRB 971214 strained theoretical models, 
GRB 990123, with a minimum redshift of 1.6 and an inferred isotropic energy release of 
1.9M0C^, takes them to the breaking point. In this paper, we present optical and infrared 
observations of GRB 990123, which both establish the distance to this luminous event, 
and through the observed break in the decay may provide the first observational evidence 
yet of beaming in a GRB. 

The optical and the infrared afterglow 



Much of our information about the parameters of GRB explosions (energy, ambient 
gas density, dynamics) has been obtained from radio and optical afterglow studies. These 
returns have motivated observers to carry out long-term multi-wavelength studies of the 
afterglow. 

We initiated an optical and IR follow-up program for GRB 990123 on 23.577 January 
1999 UT, 3.7 hr after the event, in response to the BeppoSAX preliminary localization^^. 
We used a charge coupled device (CCD) at the Palomar 60-inch telescope to image the 
field of the BeppoSAX localization. We identified a new source by comparing the first 
CCD image with an image of the field obtained from the Digital Palomar Observatory Sky 
Survey II (DPOSS), and concluded that this was the optical afterglow of GRB 990123. 
We promptly reported the discovery to the GRB Coordinates Network^^ (GCN) and the 
Central Bureau of Astronomical Telegrams. The optical source lies within the tighter 
BeppoSAX localization^^ obtained a few hours later (see Figure 1). 

Continued observations at Palomar, Keck and other observatories showed that the 
source exhibited the characteristic fading behavior seen in other GRB optical afterglows. 
The association of the optical transient (OT) with GRB 990123 is therefore secure. Table 
1 summarizes the photometric observations of the OT. Included in the table is a single 
observation^^ obtained from the Hubble Space Telescope (HST); the analysis of the HST 
data are reported elsewhere^^. Figure 2 shows the OT light curve in various bands (pri- 
marily Gunn r and K) using the photometric data shown in Table 1. 

A considerable body of literature^^'^^'^^'^^ has been developed to explain the afterglow 
phenomenon. Briefiy, afterglow can be understood as arising from synchrotron emission 
from particles shocked by the explosive debris sweeping up ambient medium. Assuming 
that the electrons behind the shock are accelerated to a power-law differential energy 
distribution with index —p, the simplest afterglow model (e.g. ref. 23) predicts that 
the afterglow fiux, f^it) oc t"z/^ where f^it) is the fiux at frequency v and t is the time 
measured with respect to the epoch of the GRB. Interpretation of the afterglow data in 
the framework of these models can, in principle, yield many interesting parameters of the 
GRB explosion (e.g. ref. 24,25). 

As the OT faded, the K band image of January 28 showed that the host galaxy was 
only a fraction of an arcsecond displaced with respect to the OT. A proper analysis of the 
decay of the OT thus requires that we correctly assess the contribution of the host galaxy to 
the fiux of the OT. The adopted values are i^host = 0.9±0.3 /uJy and rhost = 0.58±0.04 (UJy; 
see legend to Figure 2 for additional details. A discussion of the host galaxy can be found 
in Bloom et al.^^. 

In Figure 2 we present the light curve of the OT after subtracting the contribution 
from the host galaxy and correcting for Galactic extinction. No single power law can fit 
the r-band curve. We fit two power laws with a break at epoch tbreak (see Figure 2). The 
best fit parameters are air = —1.10 ± 0.03, tbreak = 2.04 ± 0.46 d and a2r = —1.65 ± 0.06. 
In contrast, the K band data are well fitted by a single power law, aK = —1.10 ± 0.11. 
We now discuss the behaviour of the afterglow in the first few days. 

The value of air and ax are similar to those measured in other afterglows (e.g. refs. 
26,27,8). In the 2-10 keV band the afterglow was observed by a variety of instruments 



aboard BeppoSAX starting 6 hr after the burst^^. From the analysis of these data, ax = 
—1.44 ± 0.07. In the framework of the afterglow models, the power law decay index 
a = 3(1 — p)/4 or a = (2 — 3p)/4, depending on whether the electrons are cooling on a 
timescale that is slower or faster than the age of the shock. A value oi p — 2.44 is 
consistent with air, «k and ax', this would require that the X-ray emitting electrons are 
in the fast cooling regime whereas the electrons emitting photons in the r and the K bands 
are in the slow cooling regime. 

As briefly summarized above, afterglow models predict that the spectrum of the OT 
should also be a power law and with the spectral index P having a specific relation to 
a. After correcting for Galactic extinction, one day after the burst we find that PrK = 
—0.8 ±0.1. This value is what is expected when the electrons emitting optical and IR 
photons are in the slow cooling regime. In contrast, comparing the X-ray flux six hours 
after the burst^^ with the interpolated r-band flux (from Figure 2) we find PrX = —0.54. 
Such a large value for f3 is inconsistent with the simple afterglow model. (Supporting this 
flat spectrum is the observation that the S-band fluxes, although limited, are comparable 
to the r-band fluxes; see Figure 2.) Selective extinction by intervening dust in the host 
galaxy is an unlikely explanation given the strong detection of the afterglow in the B band 
(Table 1) and in the U band^^. 

We draw attention to the fact that such a discrepancy between X-ray and optical 
fluxes is not uncommon and, in our opinion, a similiar discrepancy exists for GRB 971214 
and GRB 980703 as well. This topic is worthy of further study since it is indicative of 
a fundamental violation of at least one basic assumption of the simple afterglow model 
(spherical geometry, impulsive energy release and only one dominant emission mechanism 
- synchrotron emission) Indeed, later in this article, we show that the later behaviour of 
the afterglow is also not consistent with the simplest afterglow model. 

Determining the redshift 

We obtained spectra of the optical transient source with the Low Resolution Imaging 
Spectrograph^^ on the Keck II telescope on the night of January 23, 1999. These obser- 
vations consisted of three spectra, each of about 600-s duration. A quick analysis showed 
absorption features, indicating that the OT must be at or beyond a redshift of 1.61. A 
subsequent detailed data reduction yielded the spectrum displayed in Figure 3. 

The OT spectrum is marked by prominent absorption lines whose suggested identifl- 
cations are listed in Table 2, and marked in Figure 3. Similar lines are routinely seen in 
quasar spectra, and are thought to arise from absorption in intergalactic clouds. Accepting 
the identiflcations, we note that all the absorption features can be attributed to a single 
intergalactic cloud at a redshift Zabs = 1.6004 ±0.0008. The uncertainty in the redshift has 
approximately equal contributions from random and systematic errors. The absorption is 
relatively strong, suggesting the absorbing system has a high column density of gas. This 
is only the second time that an absorption redshift has been determined for a GRB OT 
(the other case being GRB 970508 [ref. 2], and possibly also GRB 980703 [ref. 9]), thus 
placing the OT unambiguously at a cosmological distance. 



Although the redshift, Zabs = 1-6, is a lower limit, it seems likely that the absorption 
originates in the GRB host galaxy. From the absence of the Lya forest in our spectrum, 
we can place a firm upper limit to the redshift of ^ < 2.9. The absence of the reported 
Lyo! forest in the spectrum"^^ obtained from the Nordic Optical Telescope (NOT)"^^ would 
further lower this limit to ^ ^ 2.2. These limits are further supported by the relatively 
strong continuum detections of the OT in the B band (see Table 1) and in the U band^^, 
which suggest little or no absorption by the intergalactic gas. 

No other convincing absorption features (aside from the normal telluric absorption) 
are detected in our spectrum in the wavelength range A ~ 4700 - 9000 A. There are no 
obvious strong emission lines out to A ~ 9600 A. The absence of other absorption systems 
in the OT spectrum is not surprising: in this redshift range there is about one metallic 
line absorber per unit redshift interval"^^. 

Our slit also included the bright galaxy ~ 10 arcsecond to the West of the OT. We 
determine a redshift of z = 0.2783 ± 0.0005 for this galaxy, in agreement with the NOT 
measurement^^. From our IR observations we note that the K band magnitude of this 
galaxy is i^ = 16.39 ± 0.03. At this redshift, this corresponds to a normal, L r^ L^ galaxy. 

The energetics of the burst 

The combination of a high fluence and a redshift, z > 1.6, imply that this burst is 
extremely energetic. For this discussion, we assume a standard Friedmann model cosmol- 
ogy with Hq = 65 km s~^ Mpc~^, Oq = 0.2, and Aq = (if Aq > then the inferred 
distance would increase, further aggravating the energetics of the burst). At the redshift 
z = 1.6004, the derived luminosity distance is D^ = 3.7 x 10^^ cm, corresponding to the 
distance modulus {m — M) = 45.39, and 1 arcsec in projection corresponds to 8.6 proper 
kpc, or 22.5 comoving kpc. 

From the observed fluence^^ (energy > 20 keV) of 5.1 x 10~^ erg cm~^ and the 
inferred D^, we estimate the 7-ray energy release, assuming the emission was isotropic, 
to be 3.4 X 10^^ erg ^ 1.9M0C^, more than the rest mass energy of a neutron star. The 
prompt optical flash detected by ROTSE"^^ had a V magnitude of 8.9, implying a peak 
luminosity in the UV band (restframe) of ~ 3.3 x 10^^ Lq. This is about a million times the 
luminosity of a normal galaxy, and about a thousand times the luminosity of the brightest 
quasars known! The prompt optical emission is attributed to emission from the reverse 
shock^'^'"^^; in contrast, the late-time afterglow which is the focus of this paper arises from 
the forward shock. 

The above energetics and the peak luminosities are truly staggering. Popular models 
suggest that GRBs are associated with stellar deaths, and not with quasars (or the nuclei 
of galaxies). Perhaps the strongest observational evidence in favor of this presumption is 
that some GRBs are found offset from their host galaxy, and are therefore not positionally 
consistent with an active nucleus. An isotropic energy release of 1.9M0C^ is, however, 
essentially incompatible with the popular stellar death models (coalescence of neutron 
stars and the currently popular model of the death of massive stars). This large energy 
release is barely consistent with exotic models such as that of baryon decay^^; in such 
models essentially the entire rest mass energy of the neutron star is released. 



There are two ways to reduce the estimated energy release. The first possibihty is that 
the burst is amphfied by gravitational lensing. Indeed, on the basis of early reports^^ of a 
possible foreground galaxy near the line of sight to the OT, we proposed that GRB 990123 
was lensed. This hypothesis received a boost when the NOT team reported^'^ possible 
intervening systems at redshifts of 0.286 and 0.201 which could be potential candidates 
for lensing. Our own subsequent, deeper imaging observations have failed to confirm any 
foreground galaxy and, as discussed above, there is little spectroscopic evidence for the 
0.286 and 0.201 absorption systems. There is therefore no direct evidence for lensing at 
this time. 

Furthermore, the lensing probability is oc A~'^ where A is the amplification provided 
by the lens. Even with A ~ 2, the expected probability for lensing at a redshift of 1.6 is 
10~^ (ref. 37). This probability is consistent with the observation that roughly 1 in 500 
cosmic radio sources, which likely have a similar redshift distribution to GRBs, are lensed. 
Thus it is not very likely^^ that one out of the 15 bursts observed by BeppoSAX so far 
would be highly magnified. 

We now consider the second possibility: the emitting surface in GRB 990123 is not 
spherical. Indeed, almost all energetic sources in astrophysics (e.g. quasars and accreting 
stellar black holes, pulsars) are not spherically emitting sources but display jet-like geome- 
try. Not surprisingly there is an extensive literature on jets in the GRB context (e.g. refs. 
39,40,41,42,43) as well as for the afterglow emission^^. 

Should the emitting surface be a jet with an opening angle of radius ^o then the 
inferred energy is reduced from the isotropic value by the beaming factor, fh ~ ^o/2- Here 
we make the reasonable assumption based on other astrophysical sources that the jet is 
two sided. Thus if 6*0 ~ 0.3 then /& ~ 5 x 10~^, and the 7-ray energy released is ~ 10^"^ 
erg, a value within reach of current models for the origin of GRBs. 

As discussed below, a marked steepening of the afterglow emission is the clearest and 
simplest signature of a jet-like geometry. The three well studied GRBs (970228, 970508, 
980703)^^'^^'^^'^^'^^'^^ exhibit, at late times (epochs longer than few days) a single power 
law decay and with indicies that are reasonable for spherical expansion, a ~ —1.2. (The 
complicated early time variations seen in afterglows of many GRBs are not relevant to the 
discussion here.) The radio afterglow of GRB 970508 has been used to argue^*^ for beaming 
in this source but this result is model dependent. It is against this backdrop that we now 
proceed to see if there is evidence for non-spherical geometry in GRB 990123. 

The shape of the burst: beaming? 

Our knowledge of the shape of the emitting region in GRBs is limited because, due 
to relativistic beaming, only a small portion (angular size ~ F"^, where F is the Lorentz 
factor of the bulk motion of the emitting material) is visible to the observer. Thus the 
observer is unable to distinguish a sphere from a jet as long as T > 9q . However, as the 
source continues its radial expansion, F will decrease, and when F < ^^^ there will be a 
marked decrease in the observed fiux. Even if the jet continues to evolve as a cone with 
a constant opening angle, the observer will see the fiux reduced by the ratio of the solid 



angle of the emitting surface (oc 9q) to that expected for a spherical expansion (oc 7 ^). 
Thus the light curve will steepen by 7^^o ^ t~^'^. This is a purely geometric effect. 

Rhoads^^ makes the important point that unless the jet is confined by some mechanism 
(as advocated in ref. 43) the jet will also expand sideways at the sound speed of shocked 
relativistic material, Cs = c/\/?>. Thus as the jet evolves, the solid angle of the emitting 
region increases as 7r(6'o + Cgtco/ct)^, where t is the time since the burst in the GRB 
frame, and tco = t/'^ is the elapsed time in the frame moving with the jet. When F 
falls below 9q the jet spreads in the lateral direction. Consequently, in this regime, F 
decreases exponentially with radius r instead of as a power law. We can therefore take 
r to be essentially constant during this spreading phase. Since t^ ~ r/2F^c, we find 

— 1/2 

F OC t^ during the spreading phase; here, Ie is the elapsed time in the frame of the 
observer at Earth. With this scaling the typical synchrotron frequency, z/^,, decreases as 
f^ and the fiux at this frequency drops as t^^. The fiux at a fixed frequency above z/^ 
drops as Fjy = Fiy^rn{^/^m)~^^~^^^'^ oc t^^, where p ~ 2.5 is the electron energy power- 
law index. This is more than one power of t steeper than that for a spherical afterglow. 

In contrast to this discussion, we expect another kind of break and that is when the 
electrons responsible for the photons in the band of interest enter the cooling regime. In 
this case, a decreases by 1/4 (e.g. ref. 23). However, unlike the break due to jet geometry 
(discussed above) this break is not broad-band. Regardless, in both cases, electron cooling 
or jets, a is expected to evolve over a timescale of tbreak- Thus the expected change in 
ct, Aa = ai — a2 is 3/4 (constant 6*0), 1 — ai/3 (spreading jet) or 1/4 (electron cooling) 
should be considered to be upper limits; see also ref. 52. 

As can be seen from Figure 2 we do indeed see a break in the r-band flux with 
Attj, = 0.55 ± 0.07; the 95% confidence interval, taking into account of covariance with 
the epoch of the break, is 0.4-0.7 (see Figure 2). On statistical grounds, we can rule out 
Aa = 1/4. The prediction of Aa = 1/4 for cooling electron cooling arises from basic 
synchrotron theory and is therefore a robust prediction. Therefore, we firmly conclude 
that the observed break is not due to electron cooling. 

The measured value of Aa support the hypothesis of a jet in GRB 990123. However, 
the lack of a break in the K-hand light curve is an issue of considerable concern. As 
explained in the legend to Figure 2, a single power-law model provides a good fit to the K- 
band data. Within statistical errors we can state that the r-band model does not conform 
to the K-band light curve. Unfortunately, the quality and quantity of the K-band data 
and the uncertainty in the adopted fiux of the host galaxy do not permit us to derive 
parameters for a broken power law model independent of those derived from the r-band 
light curve. [Future K-hand observations when the OT has essentially disappeared would 
help us resolve at least the latter uncertainty.] 

We cautiously advance the hypothesis that we are seeing evidence for a jet (perhaps 
with no sideways expansion). If so, ^0 ~ r(^break)~^ and for typical parameters 6*0 ~ 0.2 in 
which case fi, ~ 0.02. Thus the energy released in 7-rays alone is 6 x 10^^ erg. Including the 
X-ray afterglow^"^ raises the total energy release to 10^"^ erg. GRB 990123 was indeed a very 
energetic GRB. Following the discovery of the afterglow phenomenon^'^'^, it was thought 



that afterglow observations would provide only global parameters of the GRB explosion 
(e.g. energy released, ambient density), but it now appears that afterglow observations 
may give us insight into the geometry of the explosion as well. 



Acknowledgment. We are indebted to G. Neugebauer for obtaining IR data on 29 
January 1999. We thank A. Filippenko, L. Hillenbrand and J. Carpenter for kindly agreeing 
to exchange telescope time thereby enabling us to follow up this very interesting GRB and 
B. Oppenheimer for help with observations on February 7 and 8. Some of the observations 
reported here were obtained at the W. M. Keck Observatory, which is operated by the 
California Association for Research in Astronomy, a scientific partnership among California 
Institute of Technology, the University of California and the National Aeronautics and 
Space Administration. It was made possible by the generous financial support of the W. 
M. Keck Foundation. SRK's research is supported by the National Science Foundation 
and NASA. SGD acknowledges a partial support from the Bressler Foundation. 



Table 1. Photometric observations of GRB 990123 



Date 


Telescope 


Band 


Magnitude 


Reference/ 


(UT) 








Observers 


Jan 23.577 


P60 


r 


18.65 ±0.04 


SCO 


Jan 23.578 


P200 


B 


18.93 ±0.04 


LML, JSM 


Jan 23.958 


UPSO 


r 


20.00 ±0.05 


ref. 53 


Jan 24.005 


UPSO 


B 


20.16 ±0.15 


ref. 53 


Jan 24.18 


BAO 


r 


20.39 ±0.15" 


ref. 54 


Jan 24.194 


BAO 


B 


20.64 ±0.07 


ref 54 


Jan 24.547 


WO 


r 


20.93 ±0.2« 


ref 55 


Jan 24.636 


KI 


K 


18.29 ±0.04 


Koerner & Kirkpatrick 


Jan 24.934 


UPSO 


r 


21.22 ±0.08" 


ref 53 


Jan 24.978 


UPSO 


B 


22.06 ±0.20 


ref 53 


Jan 25.14 


BAO 


r 


21.37 ±0.15« 


ref 54 


Jan 25.940 


UPSO 


r 


21.73 ±0.12« 


ref 53 


Jan 26.154 


BAO 


r 


22.09 ±0.10" 


ref 54 


Jan 27.652 


KI 


K 


19.73 ±0.07 


MAM, ISM, HIT 


Jan 29.677 


KI 


K 


20.40 ± 0.08 


G. Neugebauer & LA 


Jan 30.52 


MDM 


r 


23.28 ±0.18" 


JH, lAY 


Feb 3.54 


MDM 


r 


23.88 ±0.24" 


JH, lAY 


Feb 6.6 


KI 


K 


20.84 ±0.13 


NK 


Feb 7.610 


KI 


K 


21.00 ±0.09 


SRK 


Feb 8.6 


KI 


K 


20.83 ±0.11 


SRK 


Feb 9.052 


HST 


r 


24.12 ±0.10 


ref 18 


Feb 9.6 


KI 


K 


20.97 ±0.10 


L. Hillenbrand 


Feb 10.6 


KI 


K 


21.21 ±0.11 


L. Hillenbrand 


Feb 9.654 


KII 


r 


23.91 ±0.07" 


A. Filippenko 


Feb 14.50 


MDM 


r 


24.10 ±0.10" 


JH, lAY 



Notes: 

(1) Telescope Acronyms. P60, Palomar 60-inch telescope. P200, Palomar 200-inch Hale 
telescope. UPSO, U. P. State Observatory 104-cm telescope. BAO, Bologna Astronomical 
Observatory 1.5-m telescope. WO, Fred L. Whipple Observatory 1.2-m telescope. KI, 
Keck-I 10-m telescope. MDM, 2.4-m Hiltner Telescope. KII, Keck-II 10-m telescope. 
HST, Hubble Space Telescope. 

(2) Photometric zeropoint. The night of 23 January UT was photometric at Palomar. Ab- 
solute zero-points for the Gunn-r and S-band data were obtained using the observations 
of standard star-fields^^'^^'^^. The absolute zero-point for the iiT-band observations was 
obtained using Keck-I data from the night of 24 January UT. Each bandpass zero-point 
was confirmed independently by at least two members of our group. All subsequent imag- 
ing by our group were reduced in the standard manner and photometry was propagated 
using a set of secondary stars in the GRB field. Photometry derived from data taken 
in the Cousins i?-band were recalibrated to Gunn-r; the quoted uncertainties include the 



9 



estimated zeropoint and statistical uncertainties. In this table, we include all of our pho- 
tometric data obtained to date and also data drawn from the literature which we could 
reliably place in our B, Gunn-r, and K-hand system. Magnitudes marked with the su- 
perscript a are Re magnitudes corrected to Gunn r-band using secondary reference stars 
reported in ref. 59,36. 
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Table 2. Absorption Lines Detected in the Spectrum of the OT 

Line ID ^■obs,air ^rest,vac Z Wx^obs i 

(A) (A) (A) (A) 

Al III 1862 4843.74 1862.78 1.6010 1.23 0.07 

Zn II 2026 5267.29 2026.14 1.6004 2.08 0.17 

Cr II 2062 5361.77 2062.23 1.6007 1.26 0.10 

Zn II 2062 5361.77 2062.66 1.6002 1.26 0.10 

Fell 2260 5877.17 2260.78 1.6003 0.91 0.07 

Fe II 2344 6096.14 2344.21 1.6012 3.04 0.28 

Fe II 2374 6173.87 2373.73 1.6016 3.07 0.28 

Fe II 2382 6195.29 2382.76 1.6008 3.60 0.40 

Fell 2586 6725.75 2586.64 1.6009 2.84 0.37 

Fell 2600 6759.94 2600.18 1.6005 3.48 0.14 

Mg II 2796 7269.47 2796.35 1.6003 4.59 0.30 

Mg II 2803 7289.49 2803.53 1.6008 4.80 0.52 

Mg I 2852 7416.97 2852.97 1.6005 2.87 0.18 

Notes: 

(1) Wx^obs is the observed hne equivalent width; for the restframe values, divide by (1 + ^). 

(2) The observed absorption line at 5362 A is a blend of two lines, Cr II 2062 and Zn II 
2062, and thus it appears twice. 
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Figure 1. Discovery image of the optical counterpart of GRB 990123. (Left) A 4 ar- 
cniin X 4 arcniin portion centered on the OT and extracted from the Digital Palomar 
Observatory Sky Survey II (DPOSS); the DPOSS is the digitized version of the second 
Palomar Observatory Sky Survey. (Right) A CCD image in the Gunn-r band obtained at 
the Palomar 60-inch. The new source, the presumed optical afterglow of GRB 990123, is 
marked. The circle is the 50-arcsecond localization of the X-ray afterglow^^ as obtained 
from the Narrow Field Instrument aboard BeppoSAX. Absolute astrometry on the Palo- 
mar 60-inch discovery image was obtained by comparison of 34 objects near the optical 
transient with the USNO-A2.0 V2.0 Catalogue^'^. The r.m.s. uncertainties of the astrom- 
etry are 0.28"(R.A.) and 0.26" (declination). We find the position of the optical transient 
to be, a = 15^25™30^34, 5 = +U'^45'59".ll (J2000). 
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Figure 2. Optical and infrared light curves of the transient afterglow of GRB 990123 
from 4 hours to 23 days after the burst. The observed magnitudes from Table 1 have been 
corrected for Galactic extinction and converted to flux (see below). 

(Top) The light curve of the transient + host galaxy. From ref. 61 we estimate the Galactic 
extinction in the direction of the optical transient (/, b — 73.12°, 54.64°) to be E(B— V) = 
0.01597. Thus, assuming the average Galactic extinction curve {Ry = 3.1), the extinction 
measure is Ab = 0.069, Ar = 0.041, and Ak = 0.006 mag. 

[Bottom) The inferred light curve of the transient. The contribution of the host flux 
has been subtracted from the Gunn-r and K-hand fluxes at each epoch. The adopted 
host galaxy fluxes are K = 22.1 ± 0.3 magnitude corresponding to 0.9 ± 0.3 /uJy and 
0.58 ± 0.04 |UJy (Gunn r); see below for details. The error bars include photometric error 
as well as the uncertainties in the estimated host galaxy flux. 

r-hand curve. A single power law decay model is inconsistent with the r-band data 
ix^m/d-o.f — 10). We then fit the r-band points to a broken power law model: Fj, = 
-^*(V^break)"'" for t < tbreak and F,^ = -F* (t/tbreak)"'"' for t > threak] here, t is the time in 
days since the GRB and tbreak is epoch of the break. Using a Levenberg-Marquardt x^ 
minimization^^, we find the r-band data are adequately fit by this broken power-law model 
(X^ = 12.1 for 9 d.o.f.) with the following parameters F^ = 7.0±1.9 ^uJy, tbreak = 2.04±0.46 
d, air = —1.10 ± 0.03, and a2r = —1.65 ± 0.06. The errors quoted are 1-cr confidence 
intervals for each respective parameter. Our value of air agrees very well with the Sagar 
et al.'s^^ determination (a = —1.10 ± 0.06), based on B and R observations obtained over 
the first two nights. There is some covariance between a2 and tbreak so that the confidence 
level on a subset of parameters may extend beyond the above quoted errors. For instance, 
concurrent values of (tbreak = 1 d, a2r = —1.5) are allowed at the 2-cr level as are (tbreak 
= 3 d, ct2r = —1.8). 

K-band curve. In contrast, the i^-band data can be modelled with a single power law 
index: ax = —1.12 ± 0.11 (x^ = 8.8 for 6 d.o.f.). This is an acceptable fit (probability 
of 35%). The K-hand data cannot be statistically reconciled to the best fit r-band model 
(i.e. assuming air, ct2r and tbreak to be fixed) with x^ ranging from 28 to 66 (7 d.o.f), 
depending on the value assumed for the fiux of the host galaxy. 

The fiux of the host galaxy was estimated as follows. For the r-band observations, we 
used the HST/STIS imaging where the host galaxy is clearly resolved and thus the fiux 
of the OT and the host can be measured quite accurately. In ref. 18 we discuss in great 
detail the tie between Gunn-r photometry and the STIS photometry (done with CLEAR 
filter which corresponds approximately to the V band) as well as -R-band photometry. 
The photometric ties are robust so that the fiux of the host in r band is well determined, 
0.58 ± 0.04 //Jy. 

We determined the flux of the host galaxy in K-band as follows. We used the deepest 
exposures with the best seeing images (February 9 and February 10). In these images 
the OT is clearly resolved from the host galaxy. Sets of pixels dominated by the OT or 
by the galaxy were masked, and total fluxes with such censored data were evaluated in 
photometric apertures of varying radii. Total fluxes of the OT-|-galaxy were also measured 
in the same apertures using the uncensored data. We also varied the aperture radii, and the 



13 



position and the size of the sky measurement annulus. On February 9 (February 10) UT 
we find the OT contributes 65 (57) percent (±10 percent) of the total OT+galaxy hght. 
The estimated errors of the fractional contributions of the OT to the total light reflect 
the scatter obtained from variations in the parameters of these image decompositions. In 
both epochs the fractional contribution of the host implies a flux of the host galaxy is 
0.9 ± 0.3 /uJy (i^host = 22.1 ± 0.3 mag). As a further test, we flt the K-band light curve 
assuming the flux is given as the sum of a power-law (OT) plus a constant flux (host 
galaxy). This flt yields i^host = 21.55 ± 0.20 magnitude, in agreement with the magnitude 
derived from direct imaging. However, we note that the flt is barely acceptable (x^ = 11.9, 
5 d.o.f.). 
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Figure 3. The spectrum of the optical transient of GRB 990123. Prominent absorption 
features are marked; all of these can be attributed to various absorption lines at a redshift, 
Zabs = 1.6004; the two features labeled as "atm" are telluric absorption features (A and B 
bands) of the Earth's atmosphere. 

Three 600-s integrations were taken, though the last one was terminated after 563 
s due to the onset of fog; the integrations commenced at UT 15:38, UT 15:49 and UT 
16:01, respectively. The last spectrum also suffered from a higher sky background from 
the brightening dawn sky. The sky conditions appeared to be good during the first two 
integrations, with seeing about 1 arcsecond. We used the 300 line/mm grating and the 
resulting wavelength coverage was 4700 A to 9000 A with an effective resolution of about 
11. 6A, (full width at half maximum). Spectra were taken with the 1.0 arcsecond wide 
longslit oriented at position angle (PA) 90° (i.e. East- West) and centered on the OT; the 
slit also included the prominent galaxy that is about 10 arcsecond West of the OT (see 
Figure 1). 

The usual calibration observations were not obtained given the urgency of the obser- 
vations. The calibration of the wavelength was obtained using a dispersion curve measured 
with the same grating and slit, but a somewhat different tilt, ten days earlier. Wavelengths 
of the unblended night sky emission lines were then used to adjust the zero-point of the 
wavelength solution fit; this procedure also compensated flexure of the instrument between 
the exposures. An approximate flux calibration was accomplished by using the instrument 
response curve also measured 10 days earlier; the instrument response is stable enough for 
this approach. The slit losses due to the seeing and the difference of the actual slit PA 
from the optimal (parallactic) angle at the time of observations make the net flux calibra- 
tion and the overall the shape of the spectrum somewhat uncertain. Nonetheless, we note 
that our spectroscopic flux measurement is in excellent agreement with photometry from 
interpolation of the r-band light curve. In any case, this is not important for the redshift 
measurements discussed in the text. 
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